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The distribution of times tj^N elapsed until the first j independent random walkers from a set of 
A'^ S> 1, all starting from the same site, are trapped by a quenched configuration of traps randomly 
placed on a disordered lattice is investigated. In doing so, the cumulants of the distribution of 
the territory explored by A*' independent random walkers 5jv(t) and the probability $jv(t) that no 
particle of an initial set of A*' is trapped by time t are considered. Simulation results for the two- 
dimensional incipient percolation aggregate show that the ratio between the nth cumulant and the 
nth moment of Snit) is, for large A^, (i) very large in comparison with the same ratio in Euclidean 
media, and (ii) almost constant. The first property implies that, in contrast with Euclidean media, 
approximations of order higher than the standard zeroth-order Rosenstock approximation are re- 
quired to provide a reasonable description of the trapping order statistics. Fortunately, the second 
property (which has a geometric origin) can be exploited to build these higher-order Rosenstock 
approximations. Simulation results for the two-dimensional incipient percolation aggregate confirm 
the predictions of our approach. 



I. INTRODUCTION 

Rosenstock's trapping problem is a fundamental problem of random walk theory that has been of interest for 
decades PjHiIMH- Most studies refer to the case in which a single {N = 1) random walker is placed initially at a site 
of a Euclidean or disordered lattice which is randomly filled with trap sites at a concentration c, and then performs a 
random walk until it is absorbed by one of these traps. The statistical quantity of main interest in this problem is the 
survival probability <&i(i) that the random walker is not trapped by time t, from which one can obtain the moments 
(i™) = Jg°° t™[l — ^i{t)]dt of the lifetime of this random walker. This problem has its origin in Smoluchowski's theory 
of coagulation of colloidal particles H, |^ and has been applied to many systems in physics and chemistry such 
as trapping of mobile defects in crystals with point sinks [^SQ, the kinetics of luminescent organic materials 0, 
anchoring of polymers by chemically active sites , and atomic diffusion in glasslike materials ^0] , among others. 

A generalization of the trapping problem to the case of N independent random walkers was studied by Krapivsky 
and Redner in Ref. In particular, they studied the problem of TV diffusing predators placed initially at a given 
distance from a diffusing prey in one dimension. The model in which static preys are stochastically distributed all to 
one side of the predators was the subject of a later analysis 12]. Shortly after, the order statistics of the trapping 
problem in c?-dimensional Euclidean lattices for a set of ^ 1 independent random walkers, i.e., the statistical 
description of the time fjjv elapsed until the first j out of N independent random walkers (initially starting at the 
same site) are trapped by quenched traps randomly arr ang ed on d-dimensional Euclidean lattices, was studied (and 
rigorously solved for the one-dimensional case) in Ref. |l3| . In this work, the moments (i™jv)' ~ 1,2,3, .. . were 
calculated from the probability ^j^N{t) that j random walkers of the initial set of N have been absorbed by time 
t. The key step in this calculation was the assumption of independency of the random walkers that allowed the 
establishment of a relationship between <I>j_jv(i) and the survival probability $Ar(t) = $o,Ar(0 of the full set of N 
random walkers The survival probability ^N{t) was calculated by means of Rosenstock's approximation which 
required the evaluation of the first moment of the number Si^it) of different sites visited (territory explored) by A^ 
random walkers. 

Interest in multiparticle diffusion problems has had a boost lately because of some advances in optical spectroscopy 
[l^ that make it possible to monitor events corresponding to single particles of an ensemble. The simultaneous 
tracking of iV 3> 1 fluorescently labelled particles and the analysis of the motions of the particles allows the study of 
local conditions (mechanical response, visco-elasticity) inside many complex structures such as fibrous polymers and 
the intracellular medium But biological samples and many real inorganic substances are disordered media (as 
opposed to translationally invariant Euclidean media) which are usually described as stochastic fractals 0, 0, 0| . 
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There are two main reasons for this identification: disordered systems share the statistical fractal structure of stochastic 
fractal models and diffusion is anomalous in both media. 

The single-particle {N = 1) Rosenstock trapping problem in fractal media have been thoroughly discussed by 
Blumcn, Klafter and Zumofen 2, 18]. In this paper we will study its multiparticle {N ^ 1) version, which is relevant 
for all those cases where the diffusing particles are placed (or created) in bunches. This may be especially important 
if the first or first few particles that are absorbed lead to a trigger effect. Here, we generalize to disordered fractal 
substrates the results for the order statistics of the multiparticle trapping problem obtained for Euclidean media in 
Ref. 13]. We will discover that for the two-dimensional incipient percolation aggregate, and in sharp contrast with 
the Euclidean media results, the zeroth-order Rosenstock approximation is quite incapable of describing the survival 
probabilities, and therefore the order statistics of the multiparticle trapping problem we are dealing with. This is 
because the ratio between the cumulants k„ of the distribution of S^it) and {S]y{t))™ is relatively large (and almost 
constant) for 1. We traced the origin of this behavior to the fact that the fluctuations in the number S'Ar(t) of 

distinct sites explored by a large number N of random walkers are negligible relative to the fluctuations in the number 
of sites that form the stochastic substrate. As a practical consequence, higher-order Rosenstock approximations 
are required for these media in order to give an account of the order statistics trapping problem with an accuracy 
similar to those reached by means of the zeroth-order Rosenstock approximation for Euclidean lattices. The idea 
of evaluating the survival probability for the multiparticle trapping pro blem for Euclidean media by means of the 
Rosenstock approximation was first suggested by Larralde et al. in '19'], although, to the best of our knowledge, it 
has not been implemented (except for the so-called "one-sided trapping problem" [13]) perhaps for the lack of precise 
expressions for the moments {S^{t)) of S^it)- However, for the percolation aggregate, we will discover in Sec. lHII that 
one can get a very good estimate of (5']^(t)) from the value of the first moment (S'jv(i)). This is fortunate because 
the asymptotic expansion of (5'Ar(i)) for large N is known for this medium [20]. 

The multiparticle Rosenstock trapping problem we are considering here can be seen as a stochastic generalization of 
the problem of the order statistics for the sequence of trapping times (or exit times) of a set of N independent random 
walkers, all starting form the same site at the same time, when the tra£S form a "spherical" absorbing boundary with 
a fixed radius. This problem was first studied by Lindenberg et al. ^2l\ and Weiss et al. '2^ for Euclidean lattices 
(mainly for the one dimensional case). After these pioneering works, improved results and extensions to deterministic 
and random fractal substrates have been reported [23L |2^ . A closely related multiparticle Brownian problem has 
been recently considered by Benichou et al. j25| ]: they studied the join residence times of N independent Brownian 
particles in a disc of a given radius. In particular, they studied the time spent by all N particles simultaneously in 
the disc within a given time interval, and the time which at least N — j out of N particles spend together in the disc 
within a time interval. 

The plan of the paper is as follows. In Sec. |n] the expressions that describe the order statistics of the trapping 
process are deduced. In Sec. IIIII we study the moments (S']^(i)) of the territory explored by N independent random 
walkers on two-dimensional incipient percolation aggregates by means of numerical simulation. The results of Sec. IIIII 
are applied in Sec. IIVI to obtain the survival probability $Ar(i) by means of Rosenstock's approximation. Then we 
calculate the moments {t"\) time elapsed until the first j random walkers are trapped for every j = 1,2,... 

and m = 1,2,... and compare these predictions with simulation results for the two-dimensional incipient percolation 
aggregate. A general discussion and conclusions are given in Sec. 



II. DEFINITIONS AND FUNDAMENTAL RELATIONS 



The results and definitions of this section have already been discussed in detail in the context of the trapping 
problem in Euclidean media [l^ . However, we will briefly summarize those results that are basic and necessary in 
order to follow the arguments in the rest of the paper. 

Let us first define \I'j,7v(i) as the probability that j random walkers of the initial set of N have been absorbed by 
time t by a given configuration of traps arranged on a given realization of the disordered substrate. The quantity 
of statistical interest is the average ^j^N{t) = (^j,Ar(f)) performed over all the possible outcomes of the "trapping 
experiment" carried out in a quenched configuration of traps in a given lattice realization followed by an average over 
all trap configurations and lattice realizations. We will also denote by ^'(i) the probability that a single random walker 
has not been absorbed by time t in this quenched configuration of traps placed upon a specific lattice realization. 
This is commonly known as the survival probability. It is then clear that 

vi/^.^(t) = (^^^ (1 - ^^Jy^^J^-l ^ (^^^ E(-ir^)*^"^'+" - (-i)^ (^)v^"vi/o,iv (i) 

where V^^'o,Ar(i) — X]m=o(~-'^)™(m)^o.^-™(^) J^^* backward difference formula for the jth derivative of 
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^o,N{t), d^'^o,N{t)/dNK Averaging over different configurations, and taking into account that ^o.N{t) = ^N{t) = 
(^'o,w(i)) and Ar(i) = we get 

$,-A,(i) = (-l)^(^^^V^<f>Ar(i). (2) 

Let us call Sj the state in which j particles have been absorbed and N — j particles of the initial set of N are 
still diffusing, and let hj^N (t) dt be the probability that the jth absorbed particle of the initial set of N disappears 
during the time interval {t, t + dt]. The change of the probability of Sj during the time interval (t, t + dt] is given by 
^j,N{t + dt) — $j_Ar(<). But this probability changes during this time interval by two causes: first, because the state 
Sj^i can become the state Sj if a particle of the set of iV — j + 1 particles still diffusing is trapped during the time 
interval (t,t + dt] [which happens with probability hj,N{t)dt] and, second, because Sj can become the state if 
a particle of the N — j particles still diffusing is trapped during this time interval [which happens with probability 
hj+i^Nit)dt]. Therefore ^j,N{t + dt) - $j,Ar(t) = [hj^wit) ~ hj+i^N{t)]dt, i.e., 

hj+i,N{t) = hj^N - ^ (3) 
with ft-o.A' = 0. Then, the mth moment of the time in which the jth particle is trapped is given by 

/•OC 

(t^N) = / t"^h,^^{t)dt . (4) 
Jo 

Using Eqs. igl, © and Q we find the exact recursion relation 
where 

{tT.N)=m f^-^<i>N{t)dt . (6) 







The set of Eqs. (O and ^ is remarkable because it implies that the order statistic of the trapping problem can be 
described from the knowledge of (i™Ar) only. The difference derivative in Eqs. lU and Q can be approximated 
by the ordinary derivative d' /dN^ when j <C N. This will be justified (and used) in Sec. II VI 



III. MOMENTS OF THE TERRITORY EXPLORED BY iV > 1 RANDOM WALKERS ON A 
TWO-DIMENSIONAL INCIPIENT PERCOLATION AGGREGATE 

The diffusion in percolation clusters as a model of transport in disordered media was first suggested by de Gennes 
|26j |. Percolation clusters are disordered fractals: they share the self-similarity property with deterministic fractals 
build up through deterministic rules but only in a statistical sense. In order to characterize these fractals several 
static and dynamic exponents have been defined. Perhaps, the most widely known is the fractal dimension, d/, which, 
in the case of disordered systems, is more conveniently defined using the scaling of mass with linear size, M ~ L'^f . 
However, random walkers in disordered structures are forced to follow the paths formed by the bonds between sites 
and, consequently, it is more natural to define a chemical (or topological) distance between two sites as the length of 
the shortest path along lattice bonds, If we consider the number of sites inside an hypersphere of radius V{t), 
usually known as chemical volume (also coincides with the mass if we assume that every site has a unit mass) it is 
expected that {V{ti) ~ Vb^''* , where the brackets refers to an average over all possible realization of the lattices and de, 
is the chemical dimension. Similarly, the generalized Einstein's law of diffusion for anomalous systems can be written 
in terms of the ordinary Euclidean distance, r, or the chemical one, £, and we have: 

(r2) ^ 2Dt^''^- , 

for i ^ 1, and where D and Dg are the diffusion coefficient and the chemical diffusion coefficient, respectively. The 
exponent d^ is the random walk dimension, also known as diffusion exponent. The exponent corresponding to the 
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TABLE I; Parameters appearing in the asymptotic expression of (S'jv(f)) and the ratios £2 and Ig for the two-dimensional 
incipient percolation aggregate. 







A 


c 


A 


Vo 


De 


^2 


^3 


1.65 


2.45 


1.0 


1.05 


0.8 


1.1 


0.6 


0.14 


0.015 



chemical (or topological) metric is called chemical random walk dimension. Another important exponent appearing 
in the theory of random walks in disordered media is the spectral or fracton dimension, dg — 2df /d^ = '^dg/d^ 4j. 

Rosenstock's procedure for evaluating the survival probability of a set of random walkers requires the knowledge 
of the first moments of the territory explored SN{t) by these random walkers. This is an interesting (and difficult) 
problem in itself that has already been thoroughly studied in the case of 1 independent random walkers 

tH El m HI El although only the first moment {SN{t)) has been rigourously estimated (HIHIlll. The 
average value of the territory explored by iV ^ 1 random walkers, all starting from the same site, in a disordered 
medium was analyzed in Ref. |20l and it was found that pol | 



{SN{t)) ~ SnI 



d,/2 



with 



Sn = Vo{2Di) 



di/2 



InTVy^/" 



(InTV)' 



-n \ ^ 



4") (InlnTV)^ 



(8) 



(9) 



media are known 2(1 12^ r29|. In particular, for the two-dimensional incipient percolation aggregate, Monte Carlo 
simulations in this substrate (with particles jumping from a site to one of its nearest neighbors placed at one unit 
distance in each unit time) have shown [20,|23 that di ~ 1.65, v — ~ 1) ^^^^ f^L — 2.45, ds = 2dt/d^^ ~ 1.35, 



The parameters c and Sm' are characteristic of the lattice and some of their values for several Euclidean and fractal 



-7 - In Ac'" ~ -0.62 (7 ~ 0.577216 is the Euler constant) and s\ 



(1) _ 



0.8 



Vo ~ 1.1, c ~ 1.05, 2Di ~ 1.2, 

(see Table nj. Hence we have a reasonable estimate of the asymptotic series for (5jv(i)) in Eq. © up to first order 
[n = 1), which is sufficient to account for simulation results, as Fig.^shows. 

In our simulations we also evaluated the second cumulant (variance), K2{S]s[) = {S'^) — (Sn)'^, and the third 
cumulant, ^3(6' at) = {Slj)—3{S'^){Sjsr)+2{Sj\r)^, of the territory explored 5'Ar(i) as they are necessary for implementing 
the extended Rosenstock approximation (see Sec. II V|) . We found that the ratio fc^ = Km{SN)/{SN{t))'^, although 
not very sensitive to the value of N (one can see in Figs. El and |31 that these parameters are well represented by 
k2 — 0.13 ± 0.02 and = 0.015 ± 0.02 over a wide range of N values) seems to tend to a constant value for large N 
(about 0.14 for ^2 and 0.016 for k^). This is a surprising behavior that departs considerably from that of Euclidean 
media. For example, for the d-dimensional Euclidean lattices it was found that k2 goes as 1/ln^A^ for large N. 
Also, the value of ^2 for large A^ is much smaller than for the percolation aggregate (for example, for TV = 2^°, 
k2 — 1/15^,1/30^ and 1/50^ in the one-, two- and three-dimensional Euclidean lattices, respectively E3l)i which 
has important consequences for the accuracy of Rosenstock's approximations of different orders, as we will show in 
Sec. lIVI The disordered nature of the substrate must be the reason for these remarkable differences in the behavior of 
km- What is happening is that, for large A^, the fluctuations in the number SN{t) of distinct sites explored by a large 
number A^ of random walkers are dominated by the fluctuations (over the set of stochastic lattice realizations used in 
the simulations) of the number of sites inside a hypersphere of chemical radius £ = \/2D£ t^/'^^ . We summarize this 
claim in a conjecture as follows 



lim km = hm 



(SNit)) _ Km{V) 



(10) 



where V{£) is the chemical volume (number of sites) of a hypersphere of chemical radius £ and KmiV) is the mth- 
cumulant of the distribution of V. Rigorously, the distance £ appearing in Eq. H10() is given by \fTDl t^/''™ which is the 
radius of the diffusion front in the thermodynamic limit (A^ — > 00). However, the quotient £m is not very sensitive to 
^ if a sufficiently large value of I is taken. In Fig.0|we plot a histogram for the chemical volume of a two-dimensional 
incipient percolation aggregate with £ = 100, evaluated using 2000 realizations of the lattice. Thereby we find that 
£2 — 0.14 and £3 ~ 0.015 in very good agreement with the values of ^2 and ^3, respectively, for large A^ (see figures|2 
andl^J. Consequently, we conclude that the fluctuations in Spf(t) are dominated by the disorder of the substrate 
and the influence of the value of N is completely overshadowed. Similar arguments were presented by Rammal and 
Toulouse in their pioneer work on percolation clusters |30j |. 
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FIG. 1: {SNit)) /f^"^^ versus \nN in the two-dimensional incipient percolation aggregate. The lines represent the result of the 
zeroth-order approximation (dashed line) and the first-order approximation (solid line). The symbols are simulation results 
obtained with 40000 experiments for t = 1000 (circles) and with 10000 experiments for t = 2000 (squares). 
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FIG. 2: Simulation values of k2 = K2 {Sn)/{Sn (t))^ versus In iV for the two-dimensional incipient percolation aggregate obtained 
averaging over 40000 (squares) and 20000 (circles) experiments for t = 1000. 
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FIG. 4: Histogram for the chemical volume inside a circle of chemical radius £ = 100 in the two-dimensional incipient percolation 
aggregate. A set of 2000 randomly generated clusters in a 400 x 400 square box were used in the computation. The distribution 
is clearly asymmetric around the maximum. 
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IV. ORDER STATISTICS OF THE TRAPPING PROCESS 



Assume that one has a quenched configuration of traps randomly placed on a given realization of the disordered 
lattice with probability c. If N random walkers start from an origin site free from traps at i = 0, it is clear that the 
probability that all random walkers survive by time t is given by (1 - c)^"(*). The average of this quantity over all 
possible random walks, trap configurations, and substrate realizations is known as the survival probability: $Ar(t) = 
((1 — c)"^"(*)) . Using a well-known theorem in statistics, we can define the nth-order Rosenstock approximation ^j^-* (t) 
for estimating $Ar(i) as 



<^W=cxp 



n+1 

E 



i-xy 



(11) 



where A = — ln(l — c) and kj (Sn) is the jth-cumulant of the distribution of the territory explored. In the limit n oo 



Toblem, 

mm. 



Eq. is known 
Its generalization 



we recover the exact result for $Ar(t). In the case of the single particle {N = 1) trajv 
as the extended Rosenstock approximation or truncated cumulant expansion |3, |4, 
to the N particle case was used in a one-dimensional trapping problem in Ref. 

In the previous section we showed that Monte Carlo simulations strongly suggest that K„(S'Ar) ~ /c„(S'jv(i))" for 
large N, where fc„ are constants. [Notice that ki{Sn) = (S'Ar(i)), so that fci = 1.] Therefore, inserting this result into 
Eq. Hll|) , the nth-order Rosenstock approximation becomes 



exp 



n+1 

E- 



XV 



(12) 



or equivalently, by using Eq. ©, 



$(;)(0=^exp 



ri+i 

E 



{-xy 



(13) 



We can now evaluate an approximation for the moments of the first trapping time, (^iV); by means of Eq. © 
assuming that the contribution of $Ar(t) to (i™Ar) negligible for those times for which $Ar(t) and ^'"^\t) differ 
substantially. Therefore, the substitution of Eq. H13|l into Eq. © yields 



{tT,j,)c^{tT,N)n^m i'^-^exp 
Jo 



n+1 



E 



i-xy 



1 J 



dt 



(14) 



Writing v = XSn^'''^^, the nth-order Rosenstock's estimate {t^^)n for (f^j^) becomes: 

(t™^)„ = ^(A5^)-^™/''^V„(rn) 

Us 

where 



r„(m) 



.,2m/d. 



exp 



n+1 

E 



kj 



dv 



(15) 



(16) 



Therefore, we find that the different nth-order Rosentock approximations {t^{^j^)n differ from each other only by a 
numerical factor T„(m,) (an integral) that depends only on the substrate through its spectral dimension d^ and the 
set of parameters kj, j = 1,2, . . . that come from the distribution of the chemical volume of this substrate [c.f. Eq. 
(|10|l ]. The integral in Eq. Hl()|l is trivial for n = and yields To(m) = T{2m/ds). Using the values of TableQlwe 
get ds = 2de/dlj ~ 1.375 and the estimates to(1) ~ 0.89, to(2) ~ 1.95 and to(3) ~ 10.9 for the two-dimensional 
incipient percolation aggregate. The integral in Eq. (|16|l only converges for even values of n so the next meaningful 
approximation corresponds to n = 2. Taking the values ^2 = 0.13 and k^ = 0.015 (which describe ^2 and k^ well over 
the range of values of N used in our simulations: see Figs.[21and|2Jl and evaluating the integral in Eq. (|16|l numerically, 
we found the second-order prefactors T2{m): T2(1) ~ 1.24, r2(2) w 6.04 and T2(3) w 95.6, which are systematically 
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much larger than the zeroth-order ones tq (m) , especiahy when the order m of the moment is large. This means that 
Rosenstock's approximations of order higher than zero must be necessary to provide reasonable predictions for $jv {t) 
and (i™Ar) in disordered media, especially when the moment m is large. It should be noticed that the expression for 
the first trapping time in Eq. H15() includes two approximations of different nature: (a) one due to the fact that we are 
using a finite number n of terms in the cumulant expansion [which only affects the factor T„(m)]; and (b) the other 
due to the finite number of terms considered for estimating Sn by means of the asymptotic series © . Consequently, 
it is convenient to classify these approximations by indexing them with a pair of integers (n, I): the first index gives 
the order of the Rosenstock approximation that is used, and the second gives the number of terms considered in the 
evaluation of (S'jv(t)). In this way, the approximation (n, 0) corresponds to the replacement in Eq. H15|) of Sn by the 



leading term of the series of Eq. Q , so that (t™ 



N/n 



(fi"^)„o[l + 0(l/lniV)] with 



it". 



N/nO — 



2mTn{m) 



4 [XVoi2Der^/^]"^^''^ 
Tnjm) 



IniV 



m(l~dl) 



(17) 



(18) 



and where we have absorbed all the dependence on the lattice characteristic parameters (Vq, Di, d^, . . . ) into the 
coefficient Tn(m). In the same way, if we take the two first terms in the asymptotic series of of Eq. lO, we find 
{t"\)n — (*rjv)ni[l + 0(l/lniV)^] where the approximation (n, 1) is 



N/nl 



Is' 



sJ^hniV 



IniV 



-2m/ds 



or, for In > 1, 



N/nl 



Tn{m) 



A2™M(lniV)' 



1 + midi-l)- 



,(1) 



In TV 



IniV 



(19) 



(20) 



In Fig. © we compare simulation results for the trapping time of the first particle with the theoretical predictions 
given by Eqs. (|18|) and H19I) when the parameters of Table are used. We see that the second-order Rosenstock 
approximation leads to much better results than the standard zeroth-order approximation. 

The moments j = 2, 3, . . . corresponding to the trapping of the jth particle absorbed by the traps can be 

easily estimated by means of Eq. jSJ. However, we can also obtain an explicit expression for (t™jv) if we approximate 
the difference operator V-' in ^ by the derivative /dNK The error in this approximation can be estimated from 
the equation 



VV(iV) d^f{N) 



{ANy 

In our case f{N) is (t™jv) AA^ = 1 so that 



dN^ 



d^ {tl 



O 



dN^+^ 



AN 



Nl 



dN3 



O 



Nl 



dNi^ 



(21) 



(22) 



As 



^(IniV)-- = (-l)viL-il!(ln7V)---i + }l^O[{lnN)-^-'] 



dNi 



one finds, from Eq. I^HIl [or Eq.l^HI)]: that 



dHtTN) r 



(23) 



(24) 



Taking into account that (^)/iVJ ~ 1/j! for j < N, we obtain from Eqs. jSJ, (|18|l [or H19|l ] and (|23|) the recursion 
relation 



{tT+i,N) = {tT,N) + -^{dl - l)T„(m)A-2"/<i.(inAr)-™('^l-i)-i 
3 



1 + 



1 



In AT 



(25) 
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In TV 

FIG. 5: ^(ii,Ar)A^^''°^ ™ versus In A'' for the two-dimensional incipient percolation aggregate. The lines represent the 

nth-order Rosenstock approximation that uses the /th-order approximation for 5'jv(t) with, from top to bottom, n = and 
/ — (dotted line), n = and I = 1 (dashed-dotted line), n = 2 and I = (dashed line), and n — 2 and I — 1 (solid line). In 
this and the following figures we have used fc2 — 0.13 and fea — 0.015. The symbols represent simulation results for c = 0.008 
(average over 2000 lattice realizations; circles) and c — 0.001 (20000 lattice realizations; squares). 



which can be easily solved: 



it-.) - {t\y) + m{di - 1)T„(™)A--/'^^ (iJ)ilt-i)+i 



1 + 



In TV 



where 



V'(j) = ^(i)+E- 



(26) 



(27) 



is the psi (digamma) function [s^, — —7, and 7 is the Euler constant. Equation yields 



_ T„{m) 

j,N /no - ;^2m/d,(ln7V)™(''™-l) 



by using Eq. ifT^ . and 



Tn{m) 



^ ^ ^ In TV 



(28) 



(29) 



when Eq. (|20|l is used. In Fig. © we compare the predictions for (^2.^) obtained from Eqs. (|28|l and H29|l with 
simulation results. The results are similar to that found in Fig. ||SJ) for {ti^N)- In Fig- (0 the differences {tj+i^M) — {tj.N) 
estimated from Eq. (|25l) are also plotted in a scaled form for j — 1, 2. The theoretical prediction is that, for large N , 
these points should tend to He along a straight hne (which is true) with a slope (d^ — l)T2(l)/j, i.e., a slope 0.76 for 
j = 1 and 1.0 for j = 2. The last prediction is not good for j = 2 (the fitted value is 0.90), but this should not be 
surprising because in Eq. H25|) we have ignored correction terms of order 1/lnA'^, which are very large even for huge 
values of iV. The only way to remedy this deficiency would be by increasing the number of asymptotic terms retained 



V 




FIG. 8: Simulation values of (^j,N/{tj,N) versus InA'^ for j — 1 (circles), j = 2 (squares) and j = 3 (triangles) and two trap 
concentrations: c — 0.001 (hollow symbols; 20000 lattice realizations) and c — 0.008 (filled symbols; 2000 lattice realizations) 
for the two-dimensional incipient percolation aggregate. The second order Rosenstock approximation predicts a ratio close to 
2 for large values of N. 



in the evaluation of (5jv(^)), which in turns requires knowing for n >2 [c.f. Eq. Unfortunately, these values 



of sin^ are very difficult to estimate by means of numerical simulations | 2C 
The variance of tj^N is easily obtained from Eq. (|15|l : 



and are unknown for n > 2. 



[Tn{2)-Tl{l)] 
A4/d.(ln7V)2«-i) 



1 + 2(d; 



t ^^ ^(J)+7 + 4'^ + 
f 1„ AT 



,(1) 



IniV 



In AT 



O 



\n N 



and, consequently, 



l.N 



T-„(2) 



1 



(30) 



(31) 



This is an interesting result because it means that the ratio between the variance of the first trapping time and the mean 
of that time is, for large N, independent of N. The numerical value of the ratio (Jj^n / {ti,N) is [dsr{4:/ds)/T'^{2/ds)] — 
1 ~ 1.53 for the zeroth-order Rosenstock approximation, and 2.0 for the second-order approximation. In Fig. |Slwe 
plot this ratio for j = 1,2,3 versus In N, and the second-order theoretical limit aj^^ / (ii.jv) — 2 seems to be consistent 
with the simulation data. 

Some considerations about the range of validity of the approximations developed in this paper are called for at this 
point. The approximation for {SN{t)) in Eqs. © and @ is only valid in the so called Regime II or intermediate time 
regime 20, 29]. As the integral in Eq. © was evaluated assuming that the expression for {S]y{t)) was valid for all 
times we conclude that the integral of mt™~^$Ar(t) over the short-time interval [0,tx] (^x ^ IniV being the crossover 
time between Regime I and Regime II) has to be negligible relative to (i™^), or equivalently (IniV)™ ^ (^i^at)' ^'^^ 
our approach and our results being reasonable. Taking into account the estimate for (t^jv) given in Eq. H18|l . this 
condition can be written as A <C (In A^)^''*. The concentrations of traps we have used in our simulations, c = 0.001 
and c = 0.008, verify this condition for all the values of N considered. Apart from this upper bound on c, we must 
also point out that, as also for Euclidean lattices, our results break down if most of the trapping takes place within 
the long-time Donsker-Varadhan regime ,34] . Further reference to this limitation of the theory presented in this paper 
will be made below. 
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V. SUMMARY AND CONCLUSIONS 

We have dealt with the following order statistics problem: when N independent random walkers all starting from 
the same site diffuse on a disordered lattice populated with a concentration c of static trapping sites, what is the 
distribution of the elapsed times, ij,Ar, until the first j random walkers are trapped? We were able to generalize the 
theory developed for the special case of Euclidean lattices ^3 the case of disordered substrates, and asymptotic 
expressions for the moments (iJV) with j <C iV were obtained. To this end, we used the so-called Rosenstock 
approximation, which is suitable for not very large times and small concentrations of traps, c. In this approximation 
the survival probability of the full set of N random walkers, $jv(t), is expressed in terms of the cumulants k{Sn) of 
the distribution of the territory covered 5Ar(t). 

Monte Carlo simulation results for k{Sn) in the two-dimensional incipient percolation aggregate showed that for 
large N the ratio km — Km{SN)/ {SN{t))™' with m = 2, 3 hardly depends on N and is very large in comparison with 
the corresponding Euclidean ratio. We attribute this behavior to the fluctuations in S^it) being dominated by the 
fluctuations in the volume of the medium inside a hypersphere of chemical radius £ ^ \fTDl t^l'^^ . This claim is 
supported by the result km ~ for ^ 1 found by simulations of S^ii) and V{1) in two-dimensional incipient 
percolation aggregates, where Im = K,m^)l (^(^))™ characterizes the fluctuations in the volume V(p). Therefore, the 
result km ~ Im for iV 3> 1 implies that the fluctuations in S]\[{t) are mainly accounted for by the fluctuations in V{£), 
and that the fluctuations in Sisrit) induced by the randomness of the diffusion process are irrelevant. One expects this 
also to be true for other disordered media. Hence, if {Sn{t)) is known, the cumulants of the distribution of SN{t) can 
be calculated (for any sufficiently large value of A*") from the cumulants of the distribution of the chemical volume. 
Finally, taking into account that {S]\[{t)) is reasonably well known 20], we arrive at a closed expression for the survival 
probability ^N{t) [c.f. Eq. (|12|l ] using Rosenstock's approximation. But from ^^(i) one gets the probability <!>_,■ at (t) 
that j random walkers of the initial set of N have been absorbed by time t [c.f. Eq. jSJ], so that, finally, we get the 
moments of the trapping times, (iJV) from the first moment (S'„(t)) of the territory explored! 

Comparison with simulation data shows that, in contrast with the Euclidean case, Rosenstock's approximations 
of order higher than zero are necessary to account for the order statistics results in the two-dimensional percolation 
aggregate. This is a consequence of the large value (in comparison with the Euclidean case) of km, due to the large 
fluctuations in the territory explored by the random walkers, induced, as we showed, by the spatial disorder of the 
substrate. However, some features of the order statistics of trapping hold in the disordered case: for example, we 
found that the ratio (Jj,N/{tj,N) depends only on the lattice characteristic parameters ds and km, m = 1,2, .. . for 
large N . This was conflrmed by simulations in two-dimensional percolation aggregates. 

There are some interesting problems that we still cannot answer with the theory developed in this paper. For 
example, an important quantity is the time In^n elapsed until all the particles of the inital set of N are trapped. The 
evaluation of the moments of this quantity would require specific techniques for j N as our results are limited to the 
opposite limit j ^ N. Moreover, the trapping of the last particles surely takes place in the Donsker-Varadhan time 
regime ^34j where Rosenstock's approximation for the survival probability cannot be used. The recent development 
of a Monte Carlo method to evaluate confidently the survival probability in the Donsker-Varadhan time regime for 
Euclidean lattices by Barkema et al. [s^ following a previous work of Gallos et al. Is^ could serve as starting 
point for tackling this problem. However, one should be aware that this task is not a straightforward generalization to 
disordered media of that carried out for Euclidean lattices because one has to take into consideration that, as Shapir 
[37| pointed out, the Donsker-Varandhan long-time behavior is dominated by the subset of lattice realizations that 
are more ramified (with dg = 1). Consequently, an efficient Monte Carlo technique to explore the relevant rare lattice 
realizations in percolation clusters has to be devised before. 
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